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MARKED UP SPECIFICATION 

TITLE 

A simulation and digital synthesis method of an oscillating 

phenomenon 

BACKGROUND OF THE INVENTION 
Field of th<=> Invention 

The invention refers to a digital simulation method of a 
non-linear interaction between an excitation source and a wave 
in a resonator and may be applied, in particular, to the digital 
synthesis, in real time, of an oscillating phenomenon such as 
the sound emitted by a musical instrument operating more 
particularly with sustained oscillations, such as wind or rubbed 
string instrument. 

n^nription nf the Related Art 

The phenomena of wave propagation and formation of the 
emitted sounds, in particular, by a musical instrument have 
been studied scientifically for a very long time. 

in particular, it is admitted, generally speaking, that a 
musical instrument includes, at least, one exciter, characterised 
by a non-linear characteristic, coupled possibly with certain 
linear elements (the reed, the lips, the bow, the hammer, etc..) 
and resonator elements, generally linear, where there is wave 
propagation as well as, generally, localised elements (for 
instance lateral bores or simple elements of the mass or spring 
type), generally linear as well. 

Similarly a digital instrument capable of synthes.sing the 
„ sounds emitted by a musical instrument, is composed generally 
of three main elements, respectively a first element, to sense 
the gests of a musician and to transform them into 
signals/control parameters, a second element computing the 
signal in real time, a third element converting this series of 



numbers calculated into a sound signal by means of 
digital/analogue converters, amplifier, loudspeakers. 

The present invention concerns mainly the second real 
time calculation element of the signal. 

!t is known that the digital simulation of a sound or, more 
generally, of an oscillating phenomenon, may be conducted by 
discretisation in the time domain of equations forming the 
mathematic representation of the physical phenomenon to be 
simulated. Such a model is always expressed in the form of a 
system of equations with coupled partial derivations, linear or 
non-linear. 

The simulation then consists, generally speaking, in 
computing as quickly as possible the solution of the 
acoustic/mechanical model describing the operation of the 
instrument or, at least, approximations preserving its most 
important characteristics. 

Numerous methods exist to that effect and it is possible, 
in particular, to mention modal methods (which describe the 
resonator as a resonant filter comprised of a sum of elementary 
resonances), particular methods (which describe the medium 
wherein there is a wave propagation in the form of chains of the 
type mass-springs-dampers), or the digital methods for solving 
equations with partial derivations. 

However, real time sound synthesis is difficult to realise 
and consequently, for some years, other methods have been 
developed, based on a "signal processing" formalism of the 
propagation in both directions of the resonator of the 
instrument. One may quote, for instance the methods called 
"digital wave guide" or "digital wave filter". 

Generally speaking, to represent the propagation of a 
wave, one may use, in the simplest formalism, the well known 
d'Alembert's equation, which applies to longitudinal waves 
(acoustic for instance) as well as transversal waves (vibration of 
a string for instance). In particular, in the case of the 
propagation of an acoustic wave, the acoustic pressure in all 
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points of the resonator of a wind instrument may be split into a 
si of two waves o, aooustic pressure, one propagating from 
the player to the horn, and the other from the horn to the piayer, 
which are called away-wave and return-wave, 
5 m practice, such propagation is expressed by a 

convolution equation (linear filtering), which yields the away 
wave or return wave) at one point of the resonator at each time 
I relauon to the away wave (or return wave) a, another pent at 
each mstant. In the so-called Green formulation, which may be 
,„ implemented in digital form, d'Alembert's equation spec„,e .for 
instance that this linear filter, called Green core, ,s a pure 
delay depending on the speed of propagation in the medium 

°; n 'a S sXsis mode,, these waves, respectively away and 
, 5 return waves, are represented by two signals corresponding 
respectively to both propaga.ive solutions o, the differentia, 

^""such a synthesis method is implemented, for instance in 
m . document US-A-5.332.862 which describes a synthes.se, 

comprising generally speaking: 

. one non-linear part, simulating the excter, to which two 
control parameters of the sound to be simulated are app hed^ 
These parameters are, in this case, the pressure of the player s 
breath and the pressure of his lips on the reed or the 

25 mOUthPi a e nnear part, simulating the resonator, which receives a 

. l8 „.,"„o« g„, representative c, the away w.v. .ml* r by (he 

non-linear par,, and which emits a signal noted q„ 

r»nresentative of the return wave towards said part, 
represa ntat.v ^ ^ ^ ^ ^ ^ ^ 

from the linear part and the non-linear part, 

a digital/analogue converter generating the synthes.sed 

""""obviously there are other types of synthesisers but, until 
35 now, an the methods involving a moderation (he phys.ca, 
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ohenomena within the instrument were based on the 
reposition o, the vibration inside the resonator ,n terms of 

aW ^Xi:rl;= ^ methods exhibited severe, 
shortcomings. ^ „ fQrmed , f 

instanc of everal parts of cylindrical tubes of di.feren 
diame ters the section change causes the generation of a 
diameters reflected wave at each interface. The 

transmitted wave and of a reriecieu . for 

>. . K«,^n token into account for a long time, roi 
Sr-rthrthrrar;:,' of the moderation of the voce, 

""""This type of moderation, which is identical in its 
aoorcach to the conventiona, theory of the geometrical ophcs 
r^tmpioyed, for instance. 
describe the propagation of elastic 

""""ft is Known indeed that it is interesting, in all the cases 
^en'VneTsevera, waves propagate to character se an 
■ f hv a diffusion matrix, since it is thus possmie 

interface by a amu&ivu Hiffprent waves 

access the reflexions and transmissions of the different wa 

:ZZl rrrlmur/guations are always written initiany 

, :„h physica, ^■ji:^"^ - ^ 

interface, the continuity of the pressure u, 

S,ren T, h ,rhenr often'more advantageous to use .mpedanca" 

or Stance" matrices, which „* ^^Zl 

an alternative to the scattering matrix approach J. d acoust q 

^^orthe^r'hand. when there are localised elements 

■ + h q instrument to be simulated, the 
35 other than interfaces in the instrument 10 
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■■waveguide" method must be complemented by a wave Alter 
Woe method describing such localised elements (such as mass, 
sprtnrdampers) for a correct connection between the venous 

, 5Ub " S Simi?a S rly when the acoustic resonator is composed for 
instance of a oonioa, pipe, the waves moving from the player to 
"he horn (away waves, and from .he horn to the player ^ (re urn 
waves) are different, which requires, either different 
Idolization, by means of two linear filtera correspon ,d,n, o the 
,„ Green ceres describing the propagation in each direction, or 
10 fining the cone by a success.on of shor, iength 
rviinders of different diameters. . 

Besides the sound produced by a musical instrument is 
no, derived solely from the propagation of a wave through a 
„ resonator, regardless of the complexity of its geometry, but 
esults from the non-linear coupling between said resonator and 
a„ excitation source. Such non-linear coupling ,s expressed 
physically between the physical quantities representing a cause 
pre sure in the acoustic case, strength In the meohanica ca e 
j an „ ffect (flow in the acoustic case, speed in the 
20 llZJIX called Kirchhoff variables. In the acoustic 
case this so-called Euler-Bernoulli physical law, a simplified 
version of the Navier-StoKes equations used in fluid mechanics^ 
spL ties that the acoustic pressure, at the reed or the Hps o, a 
25 'J instrument, is proportional, within one additive cohs.an, o 
the square of the acoustic flow. The formulation of the waves in 
he resonator in the form of away waves and return , waves 
• quire, therefore changing .he variables enabling .he non- 
■n coupling to be expressed, no, in relation ,o ^ Pte-ur. 
,„ flow physical variables any longer, but in relation to these new 

::; p w y ave a . ^ T - 

method. It is thus that the use of iterative or tabulation methods 
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has been recently suggested to calculate the solution of the 

non-linear system. 

Consequently, the synthesis methods used until now do 
not enable to simply express the non-linear coupling which 
exists between the excitation source and the resonator of he 
instrument and limit the physical parameterisat.on of the 
synthesis algorithms. 

pp.F F SUMMARY OP TH E INVENTION 

^7~ ntjon intends to remedy such shortcomings and to 

eliminate such limitations thanks to a new real-time simulation 
and synthesis method of an oscillating phenomenon, applicable 
es ecL.y, but without being limited thereto, to se,f ; osc,lla >ng 
wind instruments, in particular, the invention refers, 

5 simulation method enabling to take into account the physical 
process governing the operation of a real instrument and where 
the digital implementation may be particularly simple 

Moreover, from a basic method applicable to instruments 
such as a clarinet, with a cylindrical resonator, the invention 

10 m ay be adapted to the simulation of other types of wind or 

sirino instruments. 

Besides, the invention is not limited to the simulation of 
musical instruments but may be applied, generally to real-t,me 
digital synthesis of all sorts of oscillating phenomena. 
, 5 Generally, the invention therefore relates to the simulation 

*" of a non-linear interaction between an excitation source .and a 
wa ve in a resonator, by means o, a digital calculation too. from 
equations whose solution corresponds to the physical event of a 
phenomenon to be simulated. 

According to the invention, the phenomenon to be 
simulated Is translated, at each time and at a given pom. of he 
resonator by a linear relation between two variables 

t °h f e r airce-uin- 3 ;: 
35 rs= n direc h ; :r::z m 0 , . ^ « 



7 



realise a non-linear interaction between the two variables of the 
impedance or admittance relation. 

In this view, the model comprises, on the one hand, at 
least one linear part representing directly the so-called input 
s impedance or admittance of the resonator, i.e. at the point 
where the non-linear interaction occurs and, on the other hand, 
one non-linear part moderation the role of the excitat.on 
source of the phenomenon to be simulated. 

In particular, for real-time digital synthesis of an 
,o oscillating phenomenon, the invention enables, from a system of 
equations between at least two variables representative of the 
behaviour of the resonator, to establish an expression of the 
input impedance or admittance of the resonator in the form of a 
linear filter including delays, without any decomposition into 
15 away-return waves, in order to realise at least one linear part of 
the model which may be coupled with a non-linear loop involving 
the evolution of the non-linearity as expressed between the two 
variables of the impedance or admittance relat 10 n of the 
resonator. 

20 Particularly advantageously, such linear part of the model 

is composed of the sum of two elementary waveguides fulfilling 
a transfer function between the two variables of the impedance 

or admittance relation. 

According to another particularly advantageous 

25 characteristic, the model is driven by at least two parameters 
representative of the non-linear physical interaction between the 
source and the resonator, by means of a loop connecting the 
output to the input of the linear part and comprising a non-linear 
function playing the part of an excitation source for the 

30 resonator. 

Thus contrary to the synthesis methods used 
conventionally, the method according to the invention does not 
involve away and return waves, but expresses directly and 
digitally the so-called impedance, linear relation between the 
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cause and effect variables, i.e. pressure and flow m the 
acoustic case, strength and speed in the mechanical case. 

Naturally, such relation puts in evidence propagative 
elements, involving filters and delays, insofar as the physical 
phenomenon to be simulated is unchanged. 

Still thanks to the method according to the invention, 
such relation is readily processable in digital form and may then 
be associated with the non-linear relation expressed physically 
between the same variables. 

As indicated above, the present invention concerns 
therefore essentially the moderation element of a digital 
instrument which, from parameters prepared by a control 
means such as a gestural sensor operated by the player, 
computes in real time a signal liable to be transformed into a 
sound signal by a conversion element. 

In particular, for real-time synthesis, by physical 
moderation, of the sound of a musical instrument resulting 
from a non-linear coupling between the excitation source and 
the resonator, the invention enables to solve the system of 
equations representative of the phenomenon to be simulated by 
expressing directly and digitally the impedance or admittance 
linear relation between the cause and effect variables and by 
associating such linear relation in digital form with the non- 
linear relation between the same variables. Moreover, in the 
case of a resonator of a complex geometry, the former may be 
decomposed into successive elements, in order to comb.ne the 
elementary linear relations corresponding respectively to each 
element of the resonator, in order to obtain an impedance or 
admittance corresponding to the geometry of the instrument. 

As already mentioned, the invention applies, in particular, 
to real-time synthesis of the sound produced by a wind 
instrument. In such a case, the two variables of the impedance 
relation are the acoustic pressure and flow at the input of the 
resonator. 
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in the case of a cylindrical resonator having one open end, 
„ is rt cu, ar ; advantageous to realise the linear par, of he 
diattal Inscription model of the impedance equation ,n the 
, ™ o a sum of two elementary waveguides having as an 
[°Z:L souTce .he „ow a, .he input of the resonator and 
fulfilling the transfer function . ex p(-2k(«>)L) 
Ze(ro) = U5(S) = TT^iikML) 1 + exp(-2ik(co)L) 



!0 



wherein: 

. on is the wave angular frequency, 

Ze(co) is the input impedance of the resonator, 

. Pe co and Ue(.) are the Fourier transforms of the 
dimensionless values of the pressure and of the flow at 
the input of the resonator, 

1) is a function of the wave angular frequency wh.ch 
depends on the phenomenon to be simulated, 
- L is the length of the resonator 

According to another characteristic, each of the two 
waveguides involves a filter having as a transfer function: 

} _ f(co) 2 = - exp(-2ik(co)L) 

and representing a two-way trave! of a wave with the sign 
and represent! y res0 nator, each waveguide 

chanqing at the open end of the resonai 
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a ngu,ar freguency response , o, — r. ^ = — 
of the losses represented by the rmer uy 

approximated digital filter. essential characteristics 

The invention covers other essentia 

3s to waveguide models depending en the phenomenon to be 

"deed, according to an essentia, character,^ o, the 
mention, the method suggested for the s,m at, on o, ^ ™Ple 

m ore complex phenomena and, in part.cular, of d.verse types 
'"'"TnThe following desoription. we shall thus expose in detail 

- z rrrrr: 

implemented for syntheststng me clari net and then 

reed acoustic cylindrical resonator, such as ^ Z Let of 
certain adaptations for the simulation of other types 

instruments. a , r ,.k.,~o 

aEl^™ £filEIiei ^™^Ss a whole digital 
JTJ simuTaC t'Ln^ument. hy the method 

relation to the samples, with a sampling frequency 

HertZ .- * i, diaaram per a calculation by combination of 
" wav eg F :fd U ::, U^'erC fhe input impedance of a cylindrica, 

reSOn Fi 0 gure 4 provides two diagrams respectively r.pr."nt.n» 
for a cyunddea, resonator, at the top the input impedance m 
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relation to the frequency expressed in Hertz and, at the bottom, 
he angular frequency response in re.a.icn to time, in seconds 

Figure 5 is a calculation diagram of a simulate model of 
a cylindrical resonator reed-type instrument 

Fioure 6 provides two diagrams similar to Figure 3 
Figure " resonator model computed 

representing respective ly^ for a resona . , 

or^rnrdina to the invention, at tne top «hk 
"pedTnce and at the bottom the approximated impulse 

reSP0 Ftg e ure 7 provides two diagrams similar to Figure 2 
respe vely representing, for a reed model computed according 
Jo the invention, on the left the transfer function and on the nght 

thea rgrT e :ho y wr P «h n e Se variat,ons, in relation to time 
, expressed in seconds, of the internal acoustic pressure a, the 
mouthpiece of a cylindrical resonator. 

Figures 8b and 8c are enlargements of attack 

extinction transient. representi 
J Z thJ transfer function and. on the right, the angular 

according to ftte invent,^ ^ * 

epeJJof the externa, acoustic pressure respe ct.v ely^ttbe 
25 top, for a single-mode reed and, at tne oouu 

"^ng'ure 11 is a calculation diagram representing the 

-^^ie^T!r:ta;:rrzr':::= e ,e 

" ,m ^ rf ,%°r^ r Sl-. diagram representing the 
impedance of a parallel combination of cyl.ndnca, 



the 

35 resonators 
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Figure 15 provides two diagrams respectively 
representing, in the oase of a string, at the top the exact 
admittance and at the bottom the approximated admittance, m 
relation to frequency expressed in Hertz. 
5 Figure 16 is a model of a digital instrument simulating a 

strinq instrument. . 

Figure 17 shows, for a string struck, the time vanations. at 
the top, of the speed of the string at the oontaot point and, at 
the bottom, of the strength exerted by the hammer on the strmg. 

Figure 18 represents, for a string struck, the trajectory of 
the strength with time, in relation to relative displacement of the 
ham mer with respect to the string. 

Figure 19 is a general simulation diagram of an instrument 
operating by non-linear coupling between an excitabon source 
is and a resonator. 

DETAH ED D ESCRIPTION OF TH EI NVENTIQN 

iTT^TaTVe described in its application 

to a clarinet-type wind instrument. 

Figure 1 schematically represents a whole digital 
instrument for the implementation of the invention compna.ng 
" eral.y speaKing, a control element I including a gestural 
senior 1 controlled by an operator 10 and transforming the 
^tions thereof into control parameters co r , <;. Y. ' 
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actions tnereor nuu uumhv , 

moderation element U on which the control parameters act, 
Z uding pne non-linear par, 2, associated with a '-ear part 3 
and an element III creating the sound, including a mean 4 fo 
gene ating, from signals computed by the modellzat. on , , emen 
„ a signal which is transformed into sound synthes.sed by a 
^analogue converter 5. As known, the physical phenomena 
involved during the production of the sound of the Cat n. are 
expressed on the one hand, by a linear propagation equation of 
Z waves in .he pipe with loss, and on the other hand, by a 
non-linear equation linking the flow with the pressure and the 
"sp acemen? of the reed at the mouthpiece of the instrument. 
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A simulation model of the sound therefore includes a 
Hnear part of the model corresponding to the resonator ^of the 
instrument which, in the case of the clarinet is composed of a 
c nd" tube. For such geometry, assuming that the rad.u. of 
e is .arge re.ative to the thickness of the boundary 
layers, the acoustic pressure inside the tube ,s governed by an 
equation in the form: 



10 



d* 2 ' e dt z stf 

. 2 r 



-A 
i. 



a ^ R^m(^ + (cv O^J* R being the radius of the 

- ... c~ «f th P Harinet The values of the 
tube, i.e. 7mm in the case of the clarinet. . 

physical constants, in mKs units, are: c-340, lv-4.10 8, 

CP=1.4. 

5 !t=5 - 6 w°h;,e Ling for the so.utions of the exp 0(-^> 
type, co being the angular frequency of the wave, the follow.ng 
may be written: 

(1) 

k{<B ) 2 = ^-(l-iKac 2 to-/ 2 ) 

by replacing^ with the approximated value 1 2 , 
when x is small, the conventional approximated expression of 
k(co), which will be used below, becomes: 

25 k(co)- -acco 7 - 

c l 
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is written in the form of a continuous sum of all the waves likely 
to propagate in the pipe: 

p(x,t) = Jexp(-ik(<o)x)exp(ia>t)da> 

which appears as the inverse Fourier transform of the value 

^^"T^the following, the term "waveguide" will be reserved for 
the so-called Green formulation representing the propagation of 
wave in a medium, and including the dissipat.on and the 



diSPer ,n 0 this Green formalism, the transfer function of a pipe of 
length L, representing the propagation, the dissipation and the 
dispersion is : 

1 _ fa 1 . *vi - if ° L + 1 ac&T (3) 
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F(«) « expHktcJL) = exp^ac^Lj x exp^- ^+2*^2 )} 

The dissipation, represented by the modulus of Wj nd 
the dispersion represented by the phase of F(«), are therefore 
proper «oZ to I, while the propagation delay is prov.ded by 
k. The length of the pipe will be therefore the control parameter 
of the height and its radius the control parameter of the losses. 
° f We know, on the other hand, that the Fourier transform of 
the dimensionless pressures and flows at the .nput (Pe(«) 
Ue<»)> and at the output (Ps(o», U.(.)> of the resonator are 
linked by the system of equations : 

P a (o>) = cos(k((o)L)P s (co)+ / sin(k(©)L)U s (co) 
Ue(co) = /sin(k(co)L)P s (o>)+ cos(k( ( o)L)U s (co) 

Conventionally, in order to model the internal acoustic 
pressure the radiation may be neglected. The open end of the 
nstrument is therefore perfectly reflecting, which mvo.ves that 
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Ps«o)=0. This enables to express the relation between pressure 
and flow at the input of the resonator : 

P e {co) = / tan(k(co)L)U e (co) = Z e (tt)U e (<D) (4) 

where Ze(co) = i tan(k(«)L) is the normalized input 
impedance of the resonator. 

,n the case of a conventional single-mode reed or Hps 
model the dimensionless displacement x(t) of the reed w.th 
Aspect to its point of equilibrium, and the acoustic pressure 
pe(t) at the origin thereof, are linked by the equation. 
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L^^ + 5L^Ux(t) = ±PeW 



(5) 



(Or- 



2 dt 2 «> r dt 



with the sign + when the pressure tends to close the reed 
or the lips and the sign - when the pressure tends to open ,t. 

and wherein ^2«f r corresponds to the resonance 
frequency /r. for instance 2500 Hz and q r is the quality factor of 
the reed, for instance 0.2. Fniiri(3r 
By writing, the equation (5) with the s,gn + in the Fourier 
domain, we obtain the transfer function of the reed: 

Pe(ca) cOf - co 2 + iwq r co r 
whereof the angular frequency response is provided by: 



x(t) = 2^_ exp f Ico^Asin [1 ^ 



V) 
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For exemplification purposes, Figure 2 provides two 
diagram respectively indicating, on the left, the transfer 
ZZZ and on the right, the pulse response of the reed model, 
function ana, * /r= 2 500 Hz and a quality factor 

for a resonance frequency j\ *• 

q r =0,2. 
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As shown below, it should be noted that x(0)-0. 

Besides, in the case of a clarinet-type reed instrument o 
trumpet-type mouthpiece instrument, the acoustic pressure , pe(t) 
and acoustic flow ue(t) (dimensionless) at the >npu of the 
resonator are linked in a non-linear manner by the equation: 

u .(t) = \ (1 - sign(T - x(t) - 1))sign( T - p e (t)Kd - Y + x(t))^^)i ( 8 > 

in the case of a reed instrument, the parameter L is 
characteristic of the mouthpiece and takes Into . account the 
position of the Hps and the section ratio between the b.H and the 
resonator Such parameter C is proportional to the square root 
of the owning of the reed in idle position and usually ranges 

P^rlter T is the ratio between the pressure inside 
th . mouth of a player and the plating static pressure of the 
reed. For a no loss pipe, it ranges from \ for the ,„,t,.t,on of 
vibrations to \ for the position of a beating reed. 

The parameters C and y are therefore two important 
p.aying parameters insofar as they represent, respec t.vely the 
way the player pinches the reed and the pressure of the breath 

int ° "XT — in g the displacement equation of the reed 

, th. litis the impedance relation and the non-linear 
L ante Zo i appea s that the acoustic pressure and flow, a 
t mouthPieoe, are controlied hy the foliowino system of 
equations: 

J_d!xW + ^^ + x(t) = ±p e (t) W 
dt 2 © r dt 

30 n 

P (t0 ) = itanf^-V aCCl>1 ' 2LUe(£0) (10) 
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u e (t) = -0 



1(1 - sign( Y - x(t) - 1))sign( Y - PeWKO - Y + x(t))J7-P e (t)| flD 



The aim of the invention is therefore to find a formulation 
, of the impedance relation in the time domain enabling to solve 
three-equation system, by moderation the .mpedance 
relation in terms of elementary waveguides. 

To model the input impedance of the resonator m terms of 
elementary waveguides, the Fourier transform of the impedance 
10 Z e (o) is written in the form: 

sin(k(o))L) _ exp(ik(a>)L)^^pHI<W 
Z» = i tan(k(o>)L) = i ^j^T, " e xp(ik(<a)L) + exp(-ik(co)L) 

This expression may be written in the form: 

15 Pe(o) 1 (12) 

Ze(co) = UeM " T^elcpT^^kML) 1 + exp<-2ik(«,)L) 

Figure 3 represents a calculation model by combination of 
wavegu des, directly derived from such last equation and whose 
20 Z er function is the input impedance of the resonator. J 

nncoH nf a sum of two elementary waveguides. The upper 

the lower element corresponds to the second. The filter whose 

_F(©) 2 = - exp(-2ik(o>)L) re p r esents a two-way 

For exemplification purposes, Figure 4 provides two 
diagrams representing respectively, for a pipe of length L - 
0 Z and of radius R =7mm. at the top, the variation o, the input 
30 i^edance of the resonator in relation to the frequency and a 
T bottom, the angular frequency response of he 
corresponding waveguide mode,, calculafed hy inverse Four,er 
transform of the impedance. 
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The system of the three coupled physical equations (9) 
MO) ,11) enables to introduce the non-lineari.y in the form of a 
o° p nnecting .he output pe o, the resonator to ^ 
Fioure 5 provides an equivalent calculation diagram 
enaolin for the'simulation of a reed instrument or — ace 
instrument to couple in a non-linear way the displacement o 
he reed or of the hps and the acoustic pressure with the 
acoust c flow at the input o, the resonator, by computing a 
each sampled time, the interna, acoustic pressure at the 

mOU,h Thrmode, is entirely driven by the length L of the 
resonator and a. least two parameters 5 and y 

Tthe I ea par and comprising a non-linear function playing 
eTar of habitation source for the resonate. As shown on 
Fioure 4 the linear part takes up the diagram of Figure 2 and 
9 e on-ii ear function f is confrolled by both parameters . and 

input parameters, m the case or a respect 

- -- ; °- , to,h i n p eTi:y u T,siv^dr:;o;urnCsor;y 

3CC °r re^etr:; the sounds to be simulated, the 
^der^uires digital sampling and encase 
prepared, in the time domain of the q form 

to calculate the pressure pe(t at the mouthp 
the flow ue(t) but, to do so, it is necessary app 
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losses represented by the fitter F(co) by means of an 
approximated digital fitter. We shall therefore use an 

Id fictions in reiation to , h . geometry o, the resonate, in t is 
view, we express analytically the coefficients of the d.gital filter 
as functions of physical parameters. 

in practice, it is particularly advantageous to use a s ngle- 
pole filter while expressing the transfer function approximated ,n 
the form: 

~ b 0 exp(-2iTOD)^ (13) 
F(to) ~ i - a , exp(-itn) 
to 

wherein™^, /« being the sampling freguency, and 
D = f. t the pure delay corresponding to an away or return 
6 c is 

travel of the waves in the resonator. 

The parameters bo and .1 are expressed in relation to the 

|F(o>)f = |F(nr)| provided 

physical parameters so that 

20 values of a. fundamental play 

The first value adopted o1 is tnar or m= 
freguency. This enables to ensure a down-slope t,me of he 
fundamental freguency of the angular freguency response of the 
waveTutdo model using the approximated finer, identical to that 
25 of the waveguide model using the exact filter 

The second adopted value co2 is that or a n 
selected in order to obtain global rden.ica, decrease ,n the 
angular freguency responses of the waveguides, respectively, 
<=xact and approximate waveguides. 

The selection of such second value <o2 is therefore less 
30 restricted co r esponds, for Instance, to the second resonance 
pea n h case of the clarinet but. In certain cases, as shown 
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below in the case 



of the trumpet, it may be preferable to select 



a higher rank harmonic. 

The system of equations to be solved is provided by: 

|FK) 2 f (1 + a 2 - 2a i costw,)) = b 2 
|F(ffl 2 ) 2 f (1 + a 2 - 2a, cos(to 2 )) = b 2 



wherein |F(a>) 2 f = expf- 2acJ|Vj. 
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By neglecting the dispersion introduced by the non-linear 
part of the phase of F(e>), the frequencies of the harmonics are. 

- 2 i where 2k- 1 is the rank of the harmonic. 



w * = ~~Lfe 



C1 = cos( TOl ),c 2 = cosd^.F, = .F a = |F(o a ) I ,A, - F 1C1 ,A 2 

the coefficients a1 and b0 are provided by : 
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A1 - ^-Ji^Z^Lz^l^L (14) 
a, = -fT^fT 



b n = 



JaFf.fc. - c,)(A, - A a - j(Ai - A ^ - (F L l^U (15) 



J o F. - F 2 



20 From 
resonator: 



the expression of the input impedance of the 



1 _^piz?!!^M_ 

Ze(co) = lTexp(-2ik(G>)L) " 77exp(-2ik(a>)L) ' 
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and by noting z=exp(i<n), we obtain directly: 
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P e (z) 
U.(z) 



1 + 



-20 



1 + 



-2D 



„-2D 



1 - a^" 1 - b 0 z~ 
1 - a^ -1 + b 0 z 



•2D 



wherefrom the differential equation is derived: 

P e (n) = u e (n) - a lUe (n - 1) - b 0 u e (n - 2D) + a lPe (n - 1) - b oPe (n - 2D) (16) 
Figure 6 provides two diagrams respectively showing, at 
the top the variation of the input impedance of the resonator 
approximated in relation to the frequency and, at the bottom, 
the angular frequency response of the corresponding waveguide 
model, calculated from the differential equation for a cylindrical 
pipe having a length L=0.5m and a radius R=7mm. 

It can be noted that the model thus established provides a 
angular frequency response very close to that of the resonator, 
represented in Figure 2. 

As wel! as for the filter representing the losses, the 
relation between the acoustic pressure and the displacement of 
the exciter (reed or lips) must be discretized in the time domain. 
Or, the angular frequency response of the reed is a sinus wave 
damped exponentially which meets the condition x(0)=0, as 
already mentioned above. It is therefore possible to build a 
digital filter for which the displacement of the reed at the time 

t _ JL is function of the acoustic pressure at the time 
" fe 

t = n ~ 1 and not t n . This enables to verify the property x(0) = 0 
fe 

of the continuous system when the reed is subjected to a Dirac 
excitation. To comply with such condition, instead of using the 
bilinear transformation to approximate the terms ico and -co 2 , we 
use, according to the invention, the expressions ico ^y(z-z~ ) 
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and -« 2 «f e 2 (z-2 + z" 1 ), which correspond to an exact centred 

digital second order differentiation. 

With these approximations, the digital transfer function of 

the reed is provided by: 

X(z) = - 

Pe(z) & 2 + fe a {z _ 2 + z -i) + k ( Z - z~ 1 )q r co r 



CO? 2(0 r V«r J V 2 « f 
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wherefrom the differential equation is derived: 

x(n) =0x P e(n)+bi a p 6 (n-1)+a 1a x(n-1)+a 2a x(n-2) (1?) 

wherein the coefficients b1a, a1a and a2a are defined by : 

2f e 2 _ 1 fedr _ le_! 

f e 2 f e q r h 1 a2a =^-^ 

3 03 g> 2 2a>, a 0a a 0a a 0a 

Figure 7 provides, for such a reed model approximated, 
two diagrams representing, on the left the transfer function and 
on the right the angular frequency response, with a sampling 
frequency fe=44 100 Hertz, the values of the parameters being 

fr=2500Hz and qr=0,2. 

It can be seen that the diagrams obtained are very close 

to those of Figure 2 

This being established, we shall now expose an explicit 
resolution method, according to the invention, for coupling the 
differential equations with the non-linear characteristic. 

The sampled formulations of the angular frequency 
responses of the displacement of the reed and of the impedance 



23 



enable, indeed, to write the sampled equivalent of the system of 
equations (9. 10, 11) stated above, in the form: 



10 
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(18) 
(19) 
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x(n)=bi a p e (n-1 )+ai«x(n-1 )+a 2a x(n-2) 
pe(n)=u e (nH 1 u s (n-1)-boUe(n-2D)+a 1 p e (n-1)-boPe(n-2D) 

u e (n) = -<1 - sign(y - x(n) - 1))sign( Y - p e (n)Kd ~ Y + x(n))^PeH ( 20 ) 
This system of equations is implicit, because the 
calculation of p.(n) by the impedance equation requires to know 
u .(n) and such value is itself obtained from the non-linear 
equation and requires to know p.(n). Still, as already mentioned 
above the calculation of x(n) does not requires to know p.(n) 
but to know p.(n-1) which is indeed known at the time n. 

This enables, according to the invention, to solve simply 
and exactly the coupied system. To do so, the terms of the 
equations £191 and L 20l above, which do not depend on the t.me 
sample n, may, indeed, be gathered in the expressions: 

V=-a l u e (n-1)-b 0 Ue(n-2D)-f-a 1 pe(n-1)-boPe(n-2D) 
W = 1 (1 - sign( 7 - x(n) - D) * CO - Y + x(n)) 
which enables to write: 

p e (n) =u e (n)+V 

To generalise the method, it is interesting to associate u e 
(n) with a coefficient bc 0 = 1 in the case of a cylindrical pipe. 

Both equations L m and £201 above may then be written in 

the form: 



p e (n)=b 0 c 0 u e (n)+V 
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u.(n) = Wsign( 7 - PeCn^JfT 7 ^)! 



Since the term \ (1 - sign( Y - x(n) - D) cancels W when (1- 

Y+ x(n)) is negative, W always remains positive. If we 
successively consider both cases: y-pe(n)^0 and y-pe(n)<0 
corresponding respectively to the cases ue(n)>0 and ue(n)<0, 
ue(n) may be expressed exactly and without involvmg the 
unknown pe(n), in the form: 

u a (n) = - sign( Y - V)(-bc 0 W 2 + W^WfTljT^ 
Consequently, the calculation of the acoustic pressure and 
the flow at the mouthpiece at a sampled time n, may be 
conducted by using sequentially the following equations: 

x(n)=b 1a pe(n-1 )+ai a x(n-1 )+a 2a x(n-2) < 21 ) 

V=-a 1 u e (n-1)-boUe(n-2D)+a 1 p e (n-1)-boPe(n-2D) ( 22 > 

W = ~ 0 - sign(y - x(n) - - y + x{n)) < 23 > 

u (n) = 1 sign(y - V)(^bc 0 W= + W^^H^ ( 24 ) 
2 

(25) 

p e (n)=bCoU e (n)+V 

Thus the invention enables to solve in the time domain 
„ the system of equations governing the physical moderation of 
, he instrument, from a sampled formulation equivalent to the 
angular frequency response of the displacement of the reed of 
the impedance relation and of the non-linear oha ra c,er she 
which is translated into the system of equat.ons £181, £191, £201, 

30 WHere ' n |he equation (18) is a digital transcription of the model 
(m) of Figure 5, 
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. the equation (19) is a digital transcription of the 

impedance model of Figure 3, 
. the equation (20) is a digital transcription of the non- 
linear characteristic linking the displacement of the 
reed and the acoustic pressure with the acoustic flow. 
By grouping the terms which do not depend on the time 
sample n, the method according to the invention enables 
indeed, to determine the flow and the pressure at the input of 
the resonator by a sequential calculation of the equations 1211 
to (251 and to solve, in the time domain, the system of 
equations 191,(101, £1 11 governing the physical moderation of a 
clarinet-type reed instrument, in order to synthes.se the sounds 
oroduced by such an instrument. 

For exemplification purposes, Figure 7 shows the variation 
of the internal acoustic pressure at the mouthpiece, calculated 
by such a non-linear model involving waveguides, for a pipe 
having a length L=0.5m and a radius R=7mm, the values of the 
parameters being y=0A. 4=0. 4, fr=2205 Hz, qr=0. 3. 

Three phases can be observed: the attack trans ent 
corresponding to sudden increase of y and 5 . the steady state 
while T and C diminish gradually, in a linear manner, up to the 
oscillation threshold, and the extinction transient. 

In practice, the digital implementation of such a non-linear 
waveguide model, may be conducted with the use of elements 
available on the market for the gestural sensor. For in tance 
digi al implementation is possible in language C in the form of 
an externa, « clarinet » object for the environment known under 
he trade name Max-MSP, driven from MIDI controls supplied by 
a controller Yamaha WX5®. This controller measures the 
pressure of the lips on the reed, which controls the parameter s. 
and the pressure of the breath, which controls the parameter ^ 
Such information received in MIDI format (therefore between 0 
L 127) are re-standardised to correspond to the scale of he 
ph s eal parameters. The waveguide is tuned from the 
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information MIDI pitch controlled from the finger position which 
determines the length L of the pipe. 

Stiil as in a real instrument, the height changes in relation 
to the physical parameters such as y, C «>r and qr. Still, in 
reality a musical instrument is not tuned perfectly for ail the 
finger positions. The use of an all-pass filter for the 
implementation of the fractional part of the delay D is therefore 

not compulsory. 

In practice, it appears that the playing sensations of such 
, a virtual instrument are quite comparable to those of a real 
instrument. 

Stilt the method which has just been described for the 
simulation of the sound of a clarinet-type reed instrument, may 

stiil be perfected. 

s We know, indeed, that the acoustic pressure at the 

mouthpiece is not the variable representative of the sound 
perceived It is therefore interesting to calculate the external 
pressure which, for a cylindrical pipe, may be expressed as the 
time derivate of the outgoing flow: Pe *(t)=^. By neglecting 

2 o still the radiation, which involves p s (t) = 0, it becomes: 

Pe(o>) = / sin(k(co)L))U s (a>) 
U e (co) = cos(k(co)L)U s (ct>) 
wherefrom we may derive: 

25 

U s (oo) =exp(-ik(co)L)(P e (co)+U e (co)) 

From the perception viewpoint, the term exp(-ik(a>)L) is 
negligible. The expression above may therefore be simplified 
30 and becomes: 

Pe xt(t)=|(Pe(t) + Ue(t)) (26) 
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Thus from a digital viewpoint, the calculation, at each 
sampled time n, of the external pressure p 8Xt (n), is reduced to a 
simple differential between the sum of the internal pressure and 
of the flow, between the time n and the time n-1 . 

As shown in Figure 1, which represents schematically a 
whole digital instrument for the implementation of the invent.on 
in the case of a wind instrument, the signals p.(t) and u.(t) 
enabling the calculation of the external pressure p ex «<t) are 
prepared by the modelization element II from control parameters 

10 ^ For a cylindrical resonator, such modelization element II is 

of the type represented in Figure 5 and enables coupling the 
three equations (9), (10), (11)- 

The linear part 3 includes a computing block 31 of the type 
l5 represented in Figure 3, the transfer function Ze«o) of which ,s 
the input impedance of the resonator. 

The model is driven by the length L of the resonator and 
the non-linear part 2 implements a ™^*\* n *™" 
controlled by both parameters C and y and having as input 
oo parameters the pressure p.(t) calculated by the linear part 3 and 
he disp.acement x(t) of the exciter 22 calculated, ,n the case of 
the clarinet by a reed model (m) in relation to the same 
pressure p.(t) at the mouthpiece. 

From such pressure p.(t) and the flow u.(t) at the 
2 , mouthpiece, computed respectively by the linear part 3 and the 
" non-linear part 2, the bloc 4 computes the sound s.gnal p. rt (t) 
emitted by the digital instrument thanks to the converter 5. 

The method according to the invention has been described 
in detail for the simulation of a simple instrument, w.th 
cviindrical resonator, of the clarinet type. 

Still as will now be seen, the general diagram of F.gure 1 
may be applied to the simulation of more complex Phenomena. 

in particular, physical measurements have shown that the 
vibrations of a reed are more complex than a simple sine wave 
damped It has therefore appeared that the s.mple mode, 
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described above could be perfected in order to improve the 
quality of the sound produced, as perceived. To do so, we 
consider a very simplified reed model in the form of a free 
embedded string. Thus, the impedance model of a cylindrical 
pipe described above, which generates mainly odd harmonics, 
acts as a basis for the realisation of a multiple mode reed mode! 
complying with the same condition x(0) = 0. which enables to 
keep the digital resolution diagram of the equations (.211 to 1251 

specified above. 

Since the value of the impedance is real for ail the 
resonance peaks, its angular frequency response is a sum of 
cosine functions. On the other hand, it has been seen that the 
angular frequency response of a single mode reed model is a 
sine wave function. We shall therefore use, in the definition of 
the model the transfer function between a damped sine wave 
and a damped cosine wave. While keeping the same notations, 
the Fourier transform of x(t)= -^e^^}^^-^ 



is provided by 



X(cx>) = 



CO, 



cd 2 + i«q r eo r 



Similarly, the Fourier transform of 

Y (t) = , 2( ° r expf - ~ « r q r tl cosf 1 co r t is provided by 

V4^q7 V 2 ) I* J 



Y(co) = 



co r (co r q f + 2ico) 
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The transfer 

X(cp) = J±z 
Y(a>) 



function between y(t) and x(t) is then 



co r q r + 2ico 

We may thus write the transmittance model in the form: 
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X(co) = 
P e («>) 



1 - B exp 



V a J J 



1 co r q r + io | 



1 + B exp 



f ^ 



v 



(-' 



ffl a J) 
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To determine the three unknown variables ©a, C, B, three 
conditions are set. The first condition consists in keeping the 
frequency of the first peak s which is the maximum peak. To do 
so, we select (aa= car assuming that the frequency offset 
resulting from the quality factor qr can be negligible. 

The second condition, satisfied by the single mode reed 
model, consists in imposing a unit value of the transmittance 
module for © = 0 in order to maintain X(co) « Pe(co) at low 

frequencies. \ 

The third condition is an imposed value of — for the 

transmittance module at the frequency car, in order t$ r keep the 

height of the peak of the single mode reed model. 

Thanks to these three conditions, the transfer function 

thus realised reproduces the main characteristics of the single 

mode reed model. 

Both first and second conditions lead to the system of 

equations: 

2C(1 + B) J_ 



G>rA/qT+4( 1 - B ) q ' 
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2C(1 - B) = 1 
to r q r (1 + B) 

assuming: A,- * A 2 =±a) f q r and A 3 =A 2 A iqr , the 

coefficients B and C solving the system are provided by: 
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c A 3 - VA7 (28) 

WaT- "iSX 

Similarly to the case of the single mode reed, we shall 
5 prepare a digital mode! wherein the displacement of the reed at 
the sampled time t n =f is a function of the acoustic pressure, not 

at time tn but at time t^= This is possible since the angular 

frequency response of the model is the sum of damped sine 
wave functions. 

io To fulfil such condition, we use an approximation of icq in 

the form ■» - U* " D = Uexp(M - D . 

Moreover, to add to the model an additional control 
parameter, we replace the coefficient B with a filter such as 

q B b a Thus , we may adjust the damping of the 

~ 1 - a a exp(-its) 

is harmonics in relation to the damping of the fundamental. To 
keep the characteristic X(0) = 1, the parameters b a and a a are 

linked by the equation B= — — . 

The term exp[-i^] is replaced with its sampled 
equivalent: Z - Da = exp(-i^D a ) with the delay D a defined by. 
20 Da=E (^k] wherein E indicates the integer part. 



" r / 

Assuming: p = l«) r q r , the digital transfer function is 



written: 
X(z) _ 

pj& ' <P + ^(z - D)(z Da+1 - a.z°- + b a z) 

C(z~ 1 -a a z- 2 -b a z- Da ~ 1 ) 



25 = f7^(UlT^J^F^a a (P : f e )z- 2 + f.b.z-- - b a (f e - P)z-- 1 
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which leads to the differentia! equation : 

x(n)=b a1 Pa(n-1) + ba 2 Pe(n-2) + b a D lP e(n-D a -1) H-aaix{n^1) +a a2 x(n-2)H-a aD x(n-D a ) + a a oix<n-D a -1) (29) 

wherein the coefficients a.i, a a2> a aD2 , a aD i are defined by: 

and the coefficients b a i, b a 2, b a Di by: 
C . -Ca a h -Cb a 

b a i= — - , b a2 = — -— - . °aD1= 



fe' fe fe 

We notice that equation (291, thus established, enables to 
determine the dimensionless displacement x(n) of the reed at 
the sampled time n, from previous times. 

It appears therefore that, in the case of a multimode reed 
model the digital calculation diagram of the sound may be the 
same as that of the single mode reed while replacing equation 
(.211 with equation 129). which is another digital transcription of 
the model (m) of Figure 5. 

For exemplification purposes, Figure 9, which is 
analogous to Figures 2 and 7, provides, for an approximated 
multimode reed model computed according to the invention, two 
diagrams representing respectively, in solid lines, on the left, 
the transfer function and, on the right, the pulse response with a 
sampling frequency f e =44100 Hertz, the parameters having 
values f r =1837,5 Hz, q r =0,2, a a =0. 

On the same diagrams, we have superimposed, as doted 
iines, the transfer function and the angular frequency response 
of the single reed model as represented in Figure 7. 

According to another development of the invention, we 
may also improve the sound model of clarinet in order to make it 
more natural by integrating a certain noise thereto, the system 
being thus more realistic. Since the noise is created by a 
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turbulence at the reed before the beginning of the pipe, the 
noise is added to x(t). It appears on the other hand that, in 
practice, the level of noise depends on the pressure of the 
breath whereas its « colour » depends on the pressure of the 
lips on the reed. Indeed, from a physical viewpoint, the harder 
the reed is pressed, the smaller the opening between the reed 
and the pipe and the greater the turbulence. Therefore we shall 
use a simple noise model whose level is driven by y and the 

brilliance driven by 

For the digital implementation of such a noise model, we 
shall use a low-pass filter of a white noise. The transfer function 
of this filter is provided by Tb (z)-^^. the coefficient bb 

being driven by y, and the coefficient ab driven by C- The 
variation laws of bb and aa may be determined so that the 
sound simulated by the model is as realistic as possible. 

Both diagrams of Figure 10 show for exemplification 
purpose the variation of the module of the spectrum of the 
external acoustic pressure corresponding to the sound produced 
by the model, respectively on the top diagram for a single mode 
reed and on the bottom diagram for a multiple mode reed with 
additional noise, the simulation parameters being as follows: 

f r =2205 Hz, q r =0.25, a a =0, y=0.44, £=0.4, L=0,48 m, R=7.10- 3 m 

As already mentioned, the method according to the 
invention as it has just been described in details, refers to the 
simulation of sounds produced by a reed and cylindrical 
resonator musical instrument of the clarinet type. But the 
invention is not limited to such an application and may, 
conversely, be subject to numerous developments. 

Indeed from the non-linear physical model, involving 
waveguides and schematised in Figure 5, as well as of its digital 
transcription according to the sequence of equations £211 to 
£251 it will be possible, to simulate the operation of a resonator 
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having any geometry, by modifying the impedance model for 
cylindrical resonator and the associated differential equations, 
to replace the latter with impedances and more complex 
differential equations. 

While keeping the properties of the model which has just 
been described, it is possible indeed, from the general diagram 
of Figure 1, to build other complex impedance models, by 
combining certain impedance elements in parallel or in series 
and by offering digital approximations enabling explicit usage of 
the physical variables and more flexible control of the digital 
instrument. 

For exemplification purposes, we shall now describe 
certain developments of the basic model for cylindrical 
resonator, with reference to Figures 11 to 14 which represent 
equivalent calculation diagrams involving waveguides and 
corresponding to resonators having diverse geometries. 

Generally, on these diagrams which correspond to the 
computing block 31 in Figure 1, the operator C(a>) represents 
the input impedance and C~V) the input admittance of a 
cylindrical resonator, the digital model corresponding to C" 1 <co) 
being obtained by changing the sign of the coefficient b 0 only. 

A first improvement on the basic model which has just 
been described with reference to Figures 3 and 5, will enable, 
by the use of wave guides similarly, to realise a physical model 
for cylindrical resonator with terminal impedance. Such an 
element will enable, for instance, to link together parts of 
cylindrical resonators having different lengths and sections, in 
order to simulate the input impedance of a conduit of variable 
section, or still to take into account the radiation impedance. 

To do so, we consider the formalism of the transmission 
line linking the acoustic pressure and flow, respectively to the 
input of the resonator (P e (<B), U e (<o» and to its open end (P.(a>). 
Us(co)). 

Noting Z c = the impedance characteristic, we have: 
c itR 2 
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P a (o>) = cos(k(<o)L)P s {(o) + iZ c sin(k((a)L)U s ((o) 
Ue ((o) = -L sin(k(co)L)P s (G)) + cos(k(©)L)U s (co) 

By noting Z s = the output impedance and 

o/ m \ _ z c ~ Z sM we can wr jte the input impedance in two 

KW Z c +Z.(a>)' 



different ways: 



Zs(w) +itan(k(o3)L) 



P B (co) = __£c (30) 

z c u e (co) 1 + j z>2 tan(k((o)L) 



z 



c 



1 - R(g>) exp(-2ik(co)L) ^ 31) 
1 + R(co)exp(-2ik(co)L) 



The equation £311 shows therefore that the impedance of a 
cyiindricai resonator with terminal impedance may be obtained 
from the impedance of a cylindrical resonator without terminal 
is impedance, while replacing: 

exp(-2ik(co)L) with R(<o)exp(-2ik(a>)L). 

Figure 11 provides an equivalent calculation diagram 
involving waveguides, for the implementation of the equation 
{301 enabling to calculate the impedance of a cylindrical 
resonator with terminal impedance. 

Such a model enables to generate in cascade the input 
impedance of a conduit having any geometry and liable to be 
defined by a succession of elementary cylindrical conduits. 

Consequently, the invention may be applied to the 
25 simulation of the vocal conduit. 

But the invention may be subject to many other 
developments and the diagram of Figure 1 enables in particular, 
from the basic physical model for cylindrical resonator 
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schematised on Figure 5, to build specific models for the 
simulation of diverse musical instruments. 

It is thus that in a first development of the basic model, we 
shall build a model for conical resonator, usable, for instance, 
for the simulation of a saxophone. 

[f we designate by L the length of the pipe, R its input 
radius, G its opening, the distance x e between the apex of the 

cone and the input is: 
R 



x fi = 



Q 

sin — 



The input impedance, relative to the characteristic 
impedance 

tiR 2 is then provided by the expression: 



P.(m) 1 (32) 

z^t_jl_ + i 

. x_ e _ i tan(k(co)L) 
c 

which may be written in the form: 
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, , I© 

Pe(<») _ C 



1 + 



i tan(k(co)L) 

Figure 12 provides an equivalent calculation diagram 
involving waveguides, wherein the element noted D is the 
differential operator D=ico and the element noted C-1(co) 
corresponds to the diagram of Figure 2 by replacing 
-exp (-2ik(co)L) by +exp (-2ik(©)L). 

The corresponding differential equations may be 
established by following a process similar to that which has 
been described above. 
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Thus, using the bilinear transform to approximate io>, the 
digital transfer function of a conical resonator is provided by: 

P.(z) 1 w 

Z 0 U e (z) z + 1 u ^zl^lLt^sl 

"~ "" ~-1 K _-2D 

2 fe2»(z-D 1 ~ 3lZ " b ° z 
c 

By noting : Gp = 1+ —V and Gm = 1 " ^ 

2f. e 2f e 

transfer function is reduced to 



c c 



-2D-1 



P e ( Z ) 1 - (a, + 1)z- 1 + a,z- 2 - M L boz; 

Z^jT(z) = G^Ca^Gp +G m )z 1 +a 1 G m z- 2 + b 0 G m z^ 20 - b^z- 20 - 1 

wherefrom the differential equation is derived: 

p e (n) = bCoU e (n) +bc 1 u e (n-1)+bc 2 u e {n-2}+bc D Ue(n-2D)+bCDiUe(n-2D-1) 
+ ac 1 p e (n-1)+ac 2 p e (n-2)+acDP e {n-2D)+ac D iPe(n-2D-1) (33) 

wherein the coefficients bcO, bd, bc2, bcD and bcD1 are 
defined by: 



by. 



and the coefficients ad, ac2, acD and acD1 are defined 

a.,G„ + G m a,G m ^_ b 0 G m . 

ar 1_E — ac ? = -^r 21 , ac D = ■- — - — . a c oi - u q 

ac i - G p 2 G p G p 



Similarly, the invention may be applied to the case of 
short resonators which appear, for instance, in the mouthpiece 
of a brass instrument or in the bill of a reed instrument, or of a 
register hole or lateral hole. 

To do so, it is admitted that the radius of the short 
resonator is sufficient enough to keep the loss model used until 



now. 
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This approximation of a short resonator will consist in 
approximating the impedance Zl<<»)= i tan (k(co)l) for small 

values of k(m)l. 

We thus obtain the expression 

Z,(co)= i tan (k(co)l) = l(k(o>)l) = G(o>) + icoH(co) (34) 

1 - exp 0 

wherein G(co) = ) \=\ and H («>)=- (l-G(co)). 




In another development of the basic model for cylindrical 
resonator the invention also enables to simulate a more 
complex resonator, by assembling elementary impedances 
representing, in the one hand, the conduit and, on the other 
hand, the bill of a reed instrument or the mouthpiece of a brass 
instrument. 

In such a case, we shall model the mouthpiece or the bill 
by a Helmholtz resonator comprising a hemispheric cavity 
coupled with a short cylindrical pipe and a main resonator with 
conical pipe. 

The input impedance of the resonator assembly may be 
expressed by: 

. 

Ze(co) = — y- 



ka — - + 



„ 2 " . xe 

Z-j(kl(a>)M) + Z2 C ' 



1 + 



. xe 

loo — 
c 



5 



Han(k2(«)!-2) 

wherein V=|tiR^ is the volume of the hemispheric cavity, 

L1 is the length of the short pipe, L2 is the length of the conical 
pipe Z1 and Z2 are the characteristic impedances of both pipes 
which depend on their radii, k1(co) and k2(co) take into account 
the losses and of the radius R1 and R2 of each pipe. 
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This equation may be written in the form: 

x. 



103- 



iZMa>)Li) + Z 2 



ico 



1 + 



i tan(fc 2 (a>)L 2 ) 



1 . V 

1 + ICO ; 



ICO 



iZ^k^LJ + Z, 



ico- 



1 + 



itan(k 2 (a>)L 2 ) j 



We may thus establish the equivalent calculation diagram 
represented in Figure 13, wherein the operator noted Ci(co) 
corresponds to the diagram in Figure 5 and the operator noted 
S(co) corresponds to the input impedance of the conical pipe and 
to the diagram of Figure 12. 

Resolution methods similar to those which have been 
described previously, enable to express the equivalent digital 
model and the corresponding differential equations. The 
impedance of the short pipe may be modelled by using the 
approximation expressed by the previous equation (34). The 
impedance of the conical pipe is represented by the model 
corresponding to the differential equation (33). The admittance 

of the cavity ira-^r is approximated by the bilinear 

pc 2 

transformation d where d=2f e . 

z + 1 

Considering the association of the hemispheric cavity and 
of the short pipe as a Helmholtz resonator having as a 

resonance frequency co h - c , we may use such frequency o> h 

to approximate G(a>) by G(co h ) and H(<o) by H(co h ). Both 
frequencies used for the calculation of the coefficients ai and b 0 

arp m _^l^Li^!^±i^=l which correspond to the first 
are ^ - 4L(4L + 37tx e +4x e ) 



39 



impedance peak of the conical pipe, and a> 2 = co h . moreover, we 
shall use 7 = ££ = z to normalize the input impedance. 

The digital impedance of a brass-type resonator is thus 
provided by the expression: 

7 / z x - . P^_— ^ (35) 

^ z '~ Vd(z-1) 1 
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P C 2 (Z4-1) PC Ci(z) + f>C Sa(z) 

which may be simplified in the form 



Ze(Z)= ' : 



ac 0 - L k ^i ac k z ^=0 dC Dk^ 
wherefrom the differential equation may be derived : 

Pe(n)= Vbc k u e (n - k) + Xbc Dh u e (n - k - 2D) 

+ £ac k p.(n - k) + S ac DkPe(n - k - 2D) (36) 



k-1 



k-0 



The coefficients are derived from direct calculation of the 
equation (35). 

The invention may still be applied to the moderation of a 
cylindrical resonator with register holes. To that effect we use 
elements involving waveguides and corresponding respectively 
to a physical model of a cylindrical pipe with terminal 
impedance representing a pipe of length L1 between the 
mouthpiece and the register hole, a model of short pipe which 
represents the register hole of length h t and the basic model for 
cylindrical pipe representing a pipe of length L2 between the 
register hole and the open end. 
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The terminal impedance of the first part of the pipe may be 
written: 



Z s (»)= , - 



Z,C t («») Z 2 C 2 (o>) 

5 By writing such expression in the form : 

Z 2 C 2 («) 
Z t C t (co) 

we may establish a general calculation diagram by 
waveguides represented on Figure 14, which, in the case which 
l0 has ust been described, is used to calculate the term.na 
impedance by a parallel combination of the impedances Z 2 C 2 (co) 

Zt ,n1he limit case of a closed register hole, with two parts of 
the pipe having the same radius, we assume Z t Ci<«)-« 

" " d "Sing to Z s ( W )=iZ c tan(K«o)L 2 ) and Z c (.)=iZ c 
tan(k(a>)(Li+L 2 ))- 

The total input impedance of the pipe may then be 
20 expressed by : 

P (co) (C 2 (a») + C 1 (co))Z t C t (a>)j l Z^C J (©)C>^ (3?) 
Z^JJa) ~ ^lTc>)C^))Z^«(co) + Z c C 2 (co) 

Using i«=fe(1-Z- 1 ) in the impedance of the short pipe, the 
25 differential equation can be derived directly as previously. 

Thus the model for simulating the cylindrical resonator of 
the clarinet-type, obtained by direct transposition of the 
s mpi ied equations of the physical behaviour of the instrument, 
maT be adapted to the simulation of instruments with non- 
30 TyXdrL resonator, such as saxophone, trumpet or other w,nd 
instruments. 



41 



10 



15 



But the invention is not limited to such an embodiment and 
to adaptations which have just been described since, without 
departing from the scope of the claims, it may be applied to the 
simulation of other types of instruments, for instance when a 
string is rubbed as with the violin or struck as with the piano. 

Indeed, in the case of a string, it is known that using the 
formalism of mechanical transmission line, analogous to that of 
an acoustic transmission line, we establish analogous relations 
between the variables of the impedance relation. Said relations 
are the strength exerted on the mechanical element, and the 
speed of such an element resulting from such strength. Insofar 
as the quantity linked to the sound emitted is linked with the 
speed (i.e. the effect) in the mechanical case, by opposition to 
the acoustic case wherein the quantity linked with the sound 
emitted is the pressure (i.e. the cause), we prefer to describe 
the resonator in terms of admittance rather than impedance. 

The equations of the mechanical transmission line 
between a point (b) and a point (e) are then: 



V^(w) = -Jr- si rt (At (u?) L) FhioS) H cx:>H(k(^) L)VU^) 

wherein F and V represent respectively the strengths and 

speeds at each point. 

The wave number k(a>) is expressed conventionally from 
the differential equation of the movement of a string under 
J0 deflection and comprises, as in the acoustic case, propagation 
(delay), dissipation, dispersion parts (see for instance : C. 
Valette, C. Cuesta "Mecanique de la corde vibrante", Hermes, 
treatise on new technologies, serie Mecanique. 1993). 

if we assume that the end (b) of the string is fixed and the 
35 end (e) mobile, we have: 
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This relation forms the input admittance of a part of 
embedded-free string at the point where it is free, and is 
identical, within one multiplying constant, to the acoust lC 
impedance of a cylindrical resonator. It may therefore be 
represented by a diagram analogous to that of Figure 3. 

At the point where the interaction with the exciter is 
realised we write the equations of continuity between both 
parts of string 1 and 2, by considering that the total strength 
exerted is the sum of the strengths exerted on each part, while 
the speeds of each part are equal: 

/<V) = /'i M + i'M w ) 

V-'V) = V, •;..,;) = lu'U-j 

Assuming that the complete string is fixed at both ends, 
this enables to express the input admittance of the string from a 
series combination of each part of the string, i.e., in terms of 
admittance: 

v , , _ v » _ 1 = „ ~ (39) 

Such a relation is identical, within one constant, to the 
parallel combination of two cylindrical acoustic resonators and 
may therefore be represented by the diagram of Figure 14 
wherein the elements noted Ct and C2 represent the admittance 
of each part of the string. 

For exemplification purposes, Figure 15 represents, in 
relation to the frequency, at the top, the exact admittance of a 
strjng at th e eighth of its length, calculated with an expression 
of k(«) derived from a conventional model, and, at the bottom, 
the admittance approximated using an approximation of the 
losses with a digital first order filter whose coefficients are 
computed with the same method as in the acoustic case. 



43 



It should be noted that, in the case of the violin, the 
contact point between the bow and the string is very close to 
one of the ends of the string. Consequently, it is possible to use 
the « short pipe » approximation for the admittance of either of 
both parts. Moreover, insofar as the losses expressed by k(co) 
are very small in a string, it is also possible to neglect it for the 
short part. Thus, the admittance of a violin string at the contact 
point with the bow may be expressed in the same manner as the 
impedance of an acoustic conical resonator and can therefore 
be represented by the diagram of Figure 12. 

As in the case of acoustic self-oscillating instruments, the 
admittance described in this basic model comprising a string 
with two fixed ends, may be refined in order to take into account 
additional physical phenomena. The method consists again in 
associating the admittances of different elements. Thus, it is 
possible with such an approach to build an input admittance of 
a resonator composed of two strings coupled by a harmonic 
table, i.e. whose the ends are not fixed any longer, but mobile, 
indeed, in most notes of a piano, there are two or three strings 
tuned to very close frequencies, which are struck 
simultaneously by the same hammer. In such a case, the total 
admittance is expressed by an association of two admittances 
of identical strings, each of these admittances being composed 
of two parts of strings, one of these parts is expressed 
identically to the input impedance of a cylindrical pipe with 
terminal impedance. In the mechanical case, the terminal 
admittance corresponding to that of the sound board, may be 
expressed by combinations of localised elements similar to 
those employed to describe the mouthpiece or the bill (i.e. 
masses, springs, dampers), enabling to take into account one or 
several vibration modes of the sound board. 

According to the invention, the formulation of the 
resonator en terms of mechanical admittance may be 
processed, for instance in an instrument where the string is 
struck as with the piano. In such a case, and as in the acoustic 
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case, the speed of a string struck by a hammer such as that of 
a piano, may be expressed from the system with three coupled 
equations: 

/(.') = K(jj h (t) - ..;,-('«.•••'• I - ^(y k {t) - y*(t))) 

<i 2 yhit) _ i_ m (40) 



dt: 2 M'k 



- a non-linear characteristic expressing the strength in relation 
to the relative displacements and speeds of the hammer and of 
the string, 

- an equation of the dynamics of the hammer linking its 
acceleration to the strength exerted by the string thereon by 
reaction, which is analogous to the expression of the 
displacement of the reed in relation to pressure, 

- an admittance equation expressing the speed of the string in 
relation to the strength which is imposed thereupon, which is 
equivalent to the acoustic impedance relation. 

The non-linear impact characteristic used here is known 
as the Hunt-Crossley characteristic. The exponent (p) ranges 
conventionally between 2 and 3, and is not an integer. yh(n) 
designates the displacement of the hammer, ys(n) that of the 
string It should be noted that it is here a new way of writing the 
problem. Indeed, conventionally, the impedance relation used 
here is replaced with the differential equation of the movement 
of the string. 

Such a simulation method of a string instrument may be 
implemented by a digital instrument whereof the model is 
schematised on Figure 16 which is analogous to the general 
diagram of Figure 1 and wherein: 

Ye(«,) designates the input admittance of the resonator; 
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MA is, in the case of an instrument where the string is 
struck, a model of hammer, expressing its speed from the 
strength f(t); 

Vs(t) is the speed of the string and Vh(t) that of the 
hammer; 

MV is a calculation model of the speed at the bridge which 
is then radiated by the sound board, from the strength and the 
speed of the string at the hammer-string contact point; 

G is the non-linear characteristic, and groups the non- 
linear function and the calculation means of the displacements 
Yh(t) and Ys(t) from Vh(t) and Vs(t); 

Vh(0) is the control parameter acting on the block MA, 
fixing the initial speed of the hammer at impact; 

L is the control parameter of the note played. 
In the case of an instrument where the string is rubbed, G is the 
non-linear friction characteristic, whereof numerous models can 
be found in the literature and whereof the control parameters 
are the pressure of the bow on the string and its speed of 
displacement. 

In a simplified model, the block MA may be deleted. 
In the case of an instrument where the string is struck, the 
discrete time model involves, as for certain elements of the 
acoustic models, the bilinear transform to approximate the time 
derivation operators, By noting W (noted V in the acoustic case) all 
the independent terms of (n) of the differential equation linking the 
speed of the string vs(n) and the strength f(n), and vh(n) the speed 
of the hammer, the transcription of the system of equations above 
in terms of sampled signals is: 
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J(n) = K(y h (n) - |fe(«)) p 0 - <*(v h (n) - ««(«))) 
v h {n) = v h (n - 1 ) - 2^m ifin) + /( " ~ ' 



vjn) = cof(n) U 

35 I ■ •> 

y h (n) = Ph (n - 1) +• ^{vh\n) + v h {n - 1)) 

;y,(,i) = tj s (fl - 1) + ^ + - 1)) 



(41) 
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Taking into account that the exponent (p) is not an 
integer, there is no explicit solution to this system, contrary to 
the acoustic case. By substituting in the first equation the 
expressions of yh(n) and ys(n) obtained from the two last 
equations, we obtain a « fixed point » type equation: 

fir,) = (A- BfinYflC- Df(ri)) 

Although there is no analytical solution, this type of 
equation is solved conventionally using iterative methods, the 
easiest being the so-called fixed point method. The regularity in 
the time domain of the strength f(n) enables to obtain very rapid 
convergence. Indeed, we observe that two or three iterations 
are sufficient. 

The model is controlled, for the note played, by acting on 
the resonator (length, diameter, tension of the string). The 
dynamics are controlled by the initial speed of the hammer, 
obtained by fixing vh(n=0). 

For exemplification purposes of realisation, Figure 17 
represents, in relation to time, at the top the speed of the string 
at the contact point (the eighth of its length), at the bottom the 
strength exerted by the hammer on the string, solutions of the 
previous system of equations, solved by the fixed point method. 

Similarly, Figure 18 represents the time trajectory of the 
strength in relation to the relative displacement of the hammer 
with respect to the string. 

We have thus described in detail an application of the 
method according to the invention to the simulation of a string 
instrument. As for the wind instruments described previously 
and contrary to the methods known previously, the invention 
enables to dispense with away-waves and return-wave 
quantities. 

Besides, it should be noted that the simulation model of a 
string instrument, illustrated in Figure 16 is very similar to the 
wind instrument model illustrated on Figure 1. Indeed, in both 
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cases, they involve linear filters including delays, to realise 
non-iinear interaction between two physical variables, so-cailed 
Kirchhoff variables, representative of the effect and of the 
cause of the phenomenon to be simulated. 

It appears therefore that the invention may be generally 
extrapolated to the simulation of any instrument operating by 
non-linear coupling between an excitation source and a 
resonator. 

For exemplification purposes, Figure 19 represents the 
general diagram of the model of such a digital instrument 
comprising, as usual, a control element !, a modelization 
element II and an element creating the sound III. 

As previously, the modelization element ll includes a 
linear part 3 with a computing block (31) whose transfer function 
is, according to the instrument to be simulated, either the input 
impedance of the resonator Ze(co) or the admittance Ye(co) and 
one non-linear part 2 which implements a non-linear function 
21. 

The block 1 may be a gestural sensor supplying control 
parameters CL acting on the linear part 3 of the model, and 
control parameters CNL acting on the non-linear part 2. 

According to the direction of the arrows indicated on the 
block 31 on Figure 19, the linear part 3 receives from the non- 
linear part 2, from left to right, when the transfer function of the 
computing block 31 is the impedance, an effect signal E to 
produce a cause signal C which is transmitted to the non-linear 
part 2, the latter producing, from this cause signal C, a new 
effect signal E bound for the linear part 3. 

Conversely, when the transfer function of the computing 
block 31 is the admittance, the linear part 3 receives from right 
to left, from the non-linear part 2, a cause signal C and 
produces an effect signal E which is transmitted to the non- 
linear part 2 to produce a new cause signal C bound for the 
linear part 3. 
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The non-linear part 2 is associated with exciters 23 
transforming respectively the cause and effect signals to 
produce the other variables involved in the non-linear 
characteristic H. 

The block 4 includes calculation means of the sound to be 
emitted from the cause C and effect E signals, which is 
transmitted to a digital/analogue converter 5. 

The invention thus enables to simulate all sorts of 
instrument and is not limited, besides, to the field of music. 

Indeed, the method according to the invention couid be 
also applied to the simulation of other oscillating phenomena, 
thanks to an adaptation of certain differential equations and a 
selection of other non-linear characteristics and of control 
parameters taking into account the physical characteristics of 
the phenomena to be simulated. 



